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Abstract 

We study the universality issue of the transverse momentum dependent parton distributions at 
small-x, by comparing the initial/final state interaction effects in dijet-correlation in pA collisions 
with those in deep inelastic lepton nucleus scattering. We demonstrate the non-universality by 
performing an explicit calculation in a particular model where the multiple gauge boson exchange 
contributions are summed up to all orders. We comment on the implications of our results on 
the theoretical interpretation of di-hadron correlation in dA collisions in terms of the saturation 
phenomena in deep inelastic lepton nucleus scattering. 
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I. INTRODUCTION 



The Feynman parton distribution functions describe the internal structure of hadrons in 
terms of the momentum distributions of partons in the infinite momentum frame. These 
distributions depend only on the longitudinal momentum fractions of the target hadrons 
carried by the partons. The measurements of high energy hadronic processes depending 
on the Feynman parton distributions have been made possible by proving the associated 
factorization theorem, which guarantees that the parton distributions studied in different 
processes are universal [l]. 

In recent years, hadronic physicists become more interested in semi-inclusive high energy 
processes, where one can study the intrinsic transverse momentum of partons inside hadrons. 
The additional transverse momentum dependence helps to picture the parton distribution 
in a three-dimension fashion and builds the hadron tomograph through the partonic struc- 

n 

ture [2J. A number of novel hadronic physics phenomena are also strongly associated with 
the transverse momentum dependent parton distributions, such as the single transverse spin 

nil On 

asymmetries |3H6[ and small- a; saturations phenomena [7|, |8|. In the last few years, great 
progress has been made in understanding the fundamental questions associated with these 
transverse momentum dependent parton distributions, such as the gauge invariance and the 



QCD factorization 



HQ, El, 



10| . In particular, the non-universality of these distribution func- 



tions due to the final/initial state interaction effects has attracted intensive investigations. 
It has been found that the difference between final state interaction in deep inelastic scat- 
tering and the initial state interaction in Drell-Yan lepton pair production in pp collisions 
leads to an opposite sign for the single spin asymmetries in these two processes [3|, |4[. More 
complicated relation was further found for the single spin asymmetry in dijet-correlation in 



pp collisions as compared to those in DIS and Drell-Yan processes [llMlq. This eventually 
leads to a conclusion that a standard transverse momentum dependent factorization breaks 
down for this process [yj]. 

In this paper, we extend the universality discussions of the transverse momentum de- 
pendent parton distributions to the small- x domain, where the fc r dependent distributions 
have been a common practice to describe the relevant physics phenomena We expect 
the non-universality for these objects as well. However, because of differen t ap proximation 



has been made in the small- a; region, the general arguments of Refs. 14 



1 51 ] on the non- 
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universality may not apply. As far as we know, there has been no discussion on this issue 
in the literature^. The objective of this paper is to study this in detail. We will carry 
out an explicit calculation in a model where both small- X and low transverse momentum 
approximation are valid. Furthermore, we will resum the initial/final state interactions to 
all orders in perturbation to study the associated universality property. 

In particular, we investigate the universality of the small- x transverse momentum depen- 
dent parton distributions probed in hadronic dijet-correlation in nucleon-nucleus collisions, 
as compared to that in the deep elastic lepton-nucleus (nucleon) scattering. There have 
been experimental results on di-hadron correlation in dA collisions at RHIC reported by the 
STAR collaboration, and interesting phenomena were found [18j. However, the theoretical 



interpretation is not yet clear at this moment 



19|-|22|. As schematically shown in Fig. 1(a), 



two partons from the nucleon projectile and nucleus target collide with each other, and 
produce two jets in the final state, 

p + A-> Jetl + Jet2 + X , (1) 

where the transverse momenta of these two jets are similar in size but opposite to each 
other in direction. In ideal case, these two jets are produced back-to-back. However, the 
gluon radiation and intrinsic transverse momenta of the initial partons induce the imbalance 
between them. We are particularly interested in the kinematic region that the imbalance 
q± = Pi± + P 2 _i_ is much smaller than the transverse momentum of the individual jet, 
namely, \q±\ <C \Pi±\ ~ K2±|- Only in this region, the intrinsic transverse momentum 
can have significant effects. Since there are two incoming partons, both intrinsic transverse 
momenta can affect the imbalance between the two jets. For large nucleus and small-x, the 
dominant contribution should come from the intrinsic transverse momentum of the parton 
from the nucleus, for which we labeled as q± in Fig. 1(a). In the following, we will focus 
on this contributions. Of course, we emphasize that both contributions shall be taken into 
account to describe the dijet-correlation in pA scattering. 

To understand the universality property of the transverse momentum dependent parton 
distribution, we study the multi-gluon exchange between the hard scattering part with the 



1 An important factorization breaking effect has been discussed in [171 ] , which is however different from the 
non-universality issue we are investigating in this paper. See more detailed discussions below. 
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FIG. 1. (a) Schematic diagram showing that two partons from the nucleon projectile and the nu- 
cleus target collide and produce two jets in the final state, where the intrinsic transverse momentum 
q± from nucleus dominates the imbalance between the two jets; (b) illustration of initial/final state 
interactions which may affect the transverse momentum dependent quark distribution from the 
nucleus in this process; (c) as a comparison, only the final state interaction effect is present in the 
deep inelastic lepton-nucleus (nucleon) scattering. 



flfl 



nucleus target [5|, |7|, ll4j, |15|. We illustrate the generic diagrams of these interactions in 



Fig. 1(b), for the particular partonic channel qq' — > qq'. All other channels shall follow 
accordingly. Since the incoming and outgoing partons are all colored objects, there exist 
initial state interaction with the initial parton from the nucleon projectile, and final state 
interactions with the outgoing two partons. For comparison, we also plot in Fig. 1(c) the 
similar diagram for the deep inelastic lepton scattering on a nucleus target, where there is 
only final state interaction on the struck quark. Clearly, if the initial/final state interactions 
affect the transverse momentum dependence, we will conclude that they are not universal 
between these processes. We emphasize, however, that our discussions will not affect the 
universality for the transverse momentum integrated parton distributions. In particular, for 
the inclusive observables the initial/final state interaction effects can be summarized into one 
gauge link associated with the integrated parton distributions, which are universal among 
different processes [lj]. 

The rest of the paper is organized as follows. In Sec. II, we construct a model to 
investigate the universality property for the small- x parton distributions, where all gauge 
boson exchange contributions can be summed up together, including all initial and final 
state interactions. In Sec. Ill, we summarize our results, and discuss the phenomenological 
implications, in particular, on the theoretical interpretation of the di-hadron correlation in 
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dA collisions at RHIC recently observed by the STAR collaboration. 



II. INITIAL AND FINAL STATE INTERACTION EFFECTS 



We take the partonic channel qq' — > qq' as an example to show the initial/final state 
interaction effects and calculate the quark distribution in dijet correlation pA — > Jetl + 
Jet2 + X, and compare with that in the deep inelastic scattering process. At small-x, 
quark distribution is dominated by gluon splitting, and can be calculated from the relevant 

For the purpose of our calculation, we employ an Abelian model 



Feynman diagrams 



of Refs. 



15| . It is a scalar QED model with Abelian massive gluons with a mass A. 



We construct the model in such a way that the large nucleus is represented by a heavy scalar 
target with mass Ma- The scalar quarks are generated by the Abelian gluon splitting which 
is the dominant contribution at small-x. The associated quark distribution in deep inelastic 

rfl 

scattering process in this model has been calculated in [5|, |7|. 

Since we are interested in studying the final state interaction effects on the parton distri- 
bution of the nucleus, for convenience, we choose the projectile as a single scalar quark with 
charge g 2 , which differs from the charge of the scalar quark from the target nucleus Q\. In 
addition, we assume that the Abelian gluon attaches to the target nucleus with an effective 
coupling g. All the partons in this calculation is set to be scalars with a mass m. The 
coupling g-i being different from gi is to show the dependence of the parton distribution on 
the initial/final state interactions associated with the incoming parton. I 
on g 2 remains for the nucleus parton distributions, they are not universal 

We perform our calculations in the covariant gauge. The final result does not depend 
on the gauge choice. We organize the calculations in terms of orders of the coupling g. At 



the dependence 



each order, a g 
target [z,ll 



uon attaches the scalar quarks in partonic scattering part from the nucleus 
15| . As shown in Fig. [21 the lowest-order graphs contain one soft gluon exchange 
with the momentum k and k + <C P\ where the plus component of a momentum p is defined 
as p + = (p° + p z )/\/2 and the nucleus is moving in +z direction. We emphasize that 
the final result is frame-independent. We calculate the scattering amplitude in the infinite 
momentum frame of the nucleus, i.e., Pj[ — > 00. One can also perform the calculations in the 
target rest frame and take the limit of Ma — > 00, which will lead to the same result (7|. The 
small-x approximation (q + ~ k + <C P£) will be taken throughout the following calculations. 
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FIG. 2. Lowest-order graphs for di-jet production in hadron-hadron collision at small-x limit. In 
these graphs, there is one soft gluon exchange with momentum k in addition to the hard gluon 
exchange. 

We will only keep the leading contributions in this limit. In additional, we also follow the 
low transverse momentum approximation in terms of q±/Pi± (q±/Pz±) by applying the 



power counting method |13(] . An important simplification is the eikonal approximation, 
which replaces the gluon attachment to the initial and final state partons with the eikonal 
propagator and vertex. After taking the leading order contributions, we find that the q± 
dependence of these diagrams can be cast into an effective quark distribution [13], which 
takes the following form, 

= 3^/ f^^V) 2 \ Ait0t) (Ml' > ( 2 ) 
with pj_ = kj_ — qj_. Here the hard partonic part depending on hard momentum sea 
been separated from the above quark distribution in the differential cross section 
separation is only possible at the leading power contribution of q±/Pi±. The contributions 
from Fig. 2 can be written as, 

" 1 1 



e Pi i has 



13|]. This 



A^(k,p)=g gi - ' 



D x D 



(3) 



k\ + \ 2 

where we have denned D (p±) = 2xP + p~ + p 2 ± + m 2 and D\ = D(q±) and D 2 = D(p±). 
In Eq. 03]), the first term and second term in the square bracket correspond to Fig. [2] (a) 
and Fig. |2] (b), respectively. The contribution from Fig. [2] (c) and Fig. [2] (d) simply cancels. 
That means at the leading order in the coupling constant, the dependence on g 2 drops out, 
which however will change at higher orders. 

At the next-to-leading order, there are 20 graphs in total in covariant gauge. We plot 
one of these graphs as an example in Fig. ()3]) (a), and additional diagrams can be obtained 



p 



F" 



P' 



(a) 



(b) 



FIG. 3. Example diagrams for two (a) and three (b) gluons exchanges, where the gluons can attach 
all charge particles in the upper part of the diagrams from the nucleus target. 



by attaching the gluons to all incoming and outgoing scalar quarks. The total contributions 
from these diagrams are 



A® (k,p) = ^ I d[l]d[2]{gl 



1 

D~i 



1 

D~2 



1 



D 



21 



D 



22 



+ 9i92 



2 

D~2 



d 2 ki±d 2 k2A 
vJ 



kf ± +X 2 k 2 2± +\z 



D 21 

h±-k 2 ±), D u 



(4) 



D{q A 



where J d[l]d[2] stands for J r> _ _., 
k{±) and D 2 i = D(p± — ki±). Clearly, this result shows the dependence on g 2 . In order to 
check this residue dependence in the amplitude squared for the quark distribution in Eq. (2), 
we need to carry out the calculation of the amplitude up to order g 3 . 

At the g 3 order, there are 120 diagrams in total with three soft gluon-exchange (see 
e.g., Fig. 13(b)), including all possible permutation of the attachments of these three gluons 
on the target nucleus. Summing up all these graphs, we obtain the three gluon exchange 
amplitude, 



A {3) (k,p) 



3! 
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+g 2 92 



d[l]d[2]d[3] [g 3 
3 3 



1 

D~2 



3 

D~2 + D 13 



1 3 



+ 9idl 



D 21 
3 

D~2~ 



D 



21 



(5) 



D 2 i D 22 

where J d[l]d[2]d[3] follows similar definition as in Eq. (4). Again, we see the dependence 
on g 2 coupling in the second and third terms. An important cross check of these results 
is that if we set g 2 = —gi, effectively there is no charge flow in the final state, and the 
quark distribution is identical to that in the Drell-Yan process in the same model. Applying 
g 2 = —gi, we can easily see that indeed Eqs. (4,5) reproduce those calculated in Ref. [26 1. 

With the amplitude calculated up to g 3 , we are able to check the dependence on g 2 for 
the parton distribution in Eq. (2). Substituting the results in Eqs. (3-5) into Eq. (2), we find 
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that the gi dependence still remains up to order g A . If we drop all g2 terms in these results, 



a. 



we obtain the quark distribution in deep inelastic scattering process in the same model j|3, 
This clearly shows that the transverse momentum dependent quark distribution q(x, q±) is 
not universal. 

This non-universality is better illustrated when we sum up all order multi-gluon exchange 
contributions. To do that, we introduce the following Fourier transform [7|, 

A(R,r)= [ d \ e-^ R ^^A(k,p) . (6) 

K J J (27r) 2 (27r) 2 K J K J 

From the Fourier transforms of A ( - 1 ' 2 ' 3 \k,p), we can easily see that they follow the expansion 
of an exponential form, 

oo 

(R,r) = J2 A(n) ( R > r ) = iV{r L ) {l - e ^[G(i?x+rx)-G(i?x)] j e -ig 9 2G(R x ) ? ^ 

n=l 

where G(R±) = K (\R ± ) /2vr and V(r±) = K (Mr±) /2tt with M 2 = 2xP + p~ + m 2 . In 
the above result, the ^-dependence seems to only appear as a phase which may not lead to 
a physics consequence. However, because the transverse momentum q± is conjugate to the 
coordinate variable difference R±-r±, this phase will lead to a non- universality contribution 
to the quark distribution as defined in Eq. (2). Therefore, the all order result reads as, 

q( Xj q ± )=^l J dp-p- J d 2 i2 ± d 2 i2ldV ± e^<^-^) e -^ 2 ( G ^)- G ^))F(r ± )y(rl) 

x jl _ e i99i[G(R ± +r ± yG(R ± )] j jl - e^l^R'^) -G(X X )] j ^ (g) 

where r' ± = R±+r±—R' ± . This transverse momentum dependent quark distribution is clearly 
different from that calculated in deep inelastic scattering in the same model 5|, |7|]. In other 
words, this distribution is not universal. It is interesting to notice that the gi dependence 
disappears after the integration over the transverse momentum. This is consistent with the 
universality for the integrated parton distributions 



III. SUMMARY AND DISCUSSIONS 



In this paper, we have demonstrated the non-universality for the small-x parton distri- 
butions in dijet correlation, by performing an explicit scalar calculation of the initial/final 
state interaction effects, and comparing to those in deep inelastic scattering on the nucleus 



S 



target. After summing up to all orders, we find that the net effects are summarized into a 
phase which leads to a non- vanishing contribution to the quark distribution and breaks the 
universality. 

It has been argued that the light-cone gauge may simplify the factorization property for 
the hard scattering processes. For example, if we choose the advanced boundary condition 



for the gauge potential in 
state interaction effects 0, 



ight-cone gauge, the wave function of hadrons contain the final 



27|. However, as we showed in the above calculations, this does 



not help to resolve the (^-dependence in the quark distribution in the dijet correlation in 
hadronic process due to the presence of both initial and final state interactions. In other 
words, the quark distribution of the nucleus in this process has to contain the interaction 
with the incoming (outgoing) quark with charge 5% which can not be solely included into 
the wave function of the nucleus target. 

The non-universality effect found in this paper is different from the factorization breaking 
effect discussed in Ref. where the breaking effect decreases as the heavy quark mass 
(equivalent to our jet transverse momentum P\±) increases. In this paper, we are discussing 
the non-universality effects which is in the leading power, and does not vanish with large 
transverse momentum of the jet. 

Our results indicate that there might be no universality for the transverse momentum 
dependent parton distributions at small-x. This may impose a challenge to explain the dijet- 
correlation data in dA collisions at RHIC with the saturation phenomena observed in DIS 
experiments at small-x. The non-universality, on the other hand, provides an opportunity 
to study the QCD dynamics associated with the final/initial state interaction effects. These 
effects should be taken into account to understand the experimental data. We plan to 



24 



in a future publication, together with 



address this issue in a saturation model 
detailed derivation of this paper. 
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